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Abstract. To study the convergence to equilibrium in random maps we developed the spectral 
theory of the corresponding transfer (Perron-Frobenius) operators acting in a certain Banach space 
of generalized functions. The random maps under study in a sense fill the gap between expanding 
and hyperbolic systems since among their (deterministic) components there are both expanding and 
contracting ones. We prove stochastic stability of the Perron-Frobenius spectrum and developed its 
finite rank operator approximations by means of a "stochastically smoothed" Ulam approximation 
scheme. A counterexample to the original Ulam conjecture about the approximation of the SBR 
measure and the discussion of the instability of spectral approximations by means of the original 
Ulam scheme are presented as well. 

1 Introduction 

Let {Ti} be a collection of nonsingular maps from a d-dimensional smooth manifold X into itself 
with a metrics p on it, and let {pi} be a collection of nonnegative integers, such that "Y^^Pi = 1. 
A random dynamical system on the space A is a stationary stochastic process Ti, T2, • • • : A ^ A 
(i.e., with values in a space of maps), see for instance JlTt . There are several approaches in further 
detalization of this object and in the sequel we shall use the following Markov one. 

By a random map T we shall mean the Markov random process on A given by the following 
family of transition probabilities 7^(a;, A) := J^iPi'^T^'^ {A)(^) from a point a; G A to a subset A C 
A. In other words, on every time step the map Ti is chosen from the collection {Ti} independently 
from the previous choices with the probability given by the distribution {pi}. 

In the literature (especially physical) the above defined random map is often called by the 
iterated function system (with probabilities). Naturally, a pure deterministic setup was studied as 
well. This can be done as follows. According to the given collection of maps {Ti} one can define 
a new multivalued map T : A A as Tx :~ UiTiX. Assuming that the maps Ti are continuous 
and strictly contracting one can show jl^ that the multivalued map T possesses a global attractor 
Xf C A, namely the Hausdorff distance between the sets T"Y and Xf decreases exponentially 
fast for any closed nonempty set Y C X. 

It is worth note that our aim in this paper is not to study the most general setup (for example, 
one can consider an infinite (continual) collection of maps {Ti}, while the distribution {pi} might 
be place dependent (on the space variable) , but rather to give a deeper analysis of dynamical and 
statistical properties of these systems related to the convergence to equilibrium, i.e. on spectral 
problems of the corresponding transfer operators practically not studied in the literature (except 
Therefore we shall not consider well studied problems of the dimensional analysis of invariant 
sets and measures. The reader can find results and further references of this type, e.g. in j^. 

One can always realize the random map under study as a deterministic one on the extended 
phase space. Denote by O the space of one-side sequences uj :— {uji,uj2, . . .} € f2, where each uji 
belongs to the set of indices of the collection Ti. On the space f2 one defines the left-shift map 
a : ^ ^ Ct according to the rule {auj)i := Ci^i+i. The topology in the space O is defined as the 
direct product of discrete topologies, acting on the set of indices, and a Borel measure fjip on fj 
is defined as the product of distributions {pi} on the set of indices. Naturally there exists an 
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one-to-one correspondence between the sequences of maps Ti and the elements of the space Jl. On 
the extended space f2 x X one can define a new map - a skew product upon the left-shift map: 
f{u),x) := {au),T^^x). 

Let us give a few simple examples of random maps, demonstrating that even in the simplest 
cases the dynamics of random maps might be rather nontrivial. For the sake of simplicity in all our 
examples both here and in the sequel the collection of maps will consist of only two maps Ti , T2 
from the unit interval into itself and thus the corresponding distribution is described by the one 
parameter p : = pi- 

Example 1.1 Pure contractive maps: T'i(x) := x/2, T2{x) := x/2 + 1/2. 

According to ||l3| (see also the next section) for each value of the parameter < p < 1 the 
corresponding random map T has the unique invariant measure, whose support typically (when 
p ^ 1/2) is a Cantor set. 

Example 1.2 A mixed case: Ti{x) := 1 - \2x - 1|, T2{x) x/2. 

In distinction to the previous example, depending on the choice of the parameter p the properties 
of the random map T differ qualitatively. Namely for < p < ^ the random map possesses 
the unique invariant measure concentrated at 0. For p > ^ the invariant measure becomes non 
unique but there exists a unique absolutely continuous invariant one, whose density hj; is such 
that = 7fc for any k = 1,2,... and intervals Ik '■= {2~^''~'^\ 2~*'']. However the constants 7^ 

are bounded on k while 2/3<p<l/2 and go to infinity with k for < p < 2/3. 

The above examples are typical for statistical problems related to random maps and in the next 
section we shall discuss separately properties of random maps expanding on average and contracting 
on average and thus we postpone the proof of the above claims till all needed technicalities will be 
introduced. It is worth note that our results related to the expanding on average random maps are 
very similar to those known for deterministic piecewise expanding maps [Q (see also results about 
the convergence to the Sinai-Bowcn-Ruclle (SBR ) measure for multidimensional expanding on 
average random maps in [^6| ^). Therefore in the corresponding sections we mainly demonstrate 
the similarities between these two types of systems. On the other hand, spectral analysis of the 
contracting on average random maps is completely new. In fact our interest in this type of systems 
is due to the analysis of spectral properties of Anosov maps in Q, where the presence of stable 
foliations relates the situation to our case. Some of the methods and ideas used in this paper 
were originated from the construction in Q needed to study the behavior of the transfer operators 
'along the stable foliation' and are based technically on the establishing of the so called Lasota- 
Yorke type inequalities. The nature of the system under study gives some advantages compare to 
the Anosov maps, in particular, the functional Banach spaces in our case do not depend on the 
finite structure of the map, the spectrum stability results are proven for a much broader class of 
random perturbations, and a more direct approximation scheme is elaborated. 

The paper is organized as follows. In the following two sections (2 and 3) we introduce the 
notions of contractive and expanding on average random maps and discuss some of their basic 
properties slightly generalizing known results about them. The main results of the paper are 
described in section 4, where we construct and study the Perron-Frobenius spectrum for random 
maps. We prove also stochastic stability of isolated eigenvalues in this spectrum and construct by 
means of a special random smoothing two schemes of their finite rank operator approximation. To 
some extent these schemes generalize the well known Ulam approximation scheme (see also 
discussion of its realization in jl| |^, |ll], |l8|). We discuss in detail why the spectrum approximation 
by means of the original Ulam scheme is not stable in the case of random maps (see also |^). 
Since the original Ulam conjecture about the approximation of the SBR measure (the leading 



element of the spectrum) still holds in our setting, in Lemma 4.13 we construct the first (to the 
best of our knowledge) counterexample to this conjecture. Note however that the map in this 
counterexample is not only discontinuous, but the discontinuity occurs in a periodic turning point 
(compare to instability results for general random perturbations in Q). 



2 Contracting on average random maps 



An important and well known case of random maps corresponds to the situation when all the maps 
Ti are continuous and their contracting constants 

\T{x)-T{y)\ . 
At- := sup J— < A < 1. 

x,yex P(x,yj 

Only the list of literature dedicated to various theoretical and applied aspects of the theory of this 
class of random maps would fill several pages, therefore we only give a reference to one of the first 
(and giving a very essential contribution) works in this field - jl^, to the monograph and to 
one of the recent publications where the reader can find more references to recent results. 
Note that these papers are dedicated mainly to questions related to dimensional and multifractal 
properties of invariant sets and measures which we shall not touch in this work. 

We start the analysis of ergodic properties of these systems using a weaker assumption, namely 
that the random map T is contracting on average. 

Definition 2.1 We shall say that the random map T is contracting on average if its contracting 
constant 

A^ti^^P.At, < 1. 

i 

Let A4 be the space of probability measures on X. Then the Markov operator associated with 
the random map T for each probabilistic measure fi E M can be written as 

TtJ. ■■= ^PifioT^-^. 

i 

Following the standard scenario we introduce in the space M the special metrics (often called 
Hutchinson metrics) JlSf : 

p/f(/i, i^) := sup hdii — J hdiy, h G C'^(A,Il), \h{x) — h{y)\ < p{x — y), ^x,y G , 

where p is the metrics on our original phase space X. It is well known (see, for example, that 
the pair [AA^pfj) defines a compact metric space. 

Lemma 2.1 Let h : X ^ he a continuous Junction and let p, E A4. Then 

J hd{Tp)=Y,Pr J hoT.dp. 

Proof. For each continuous function h : X H there is a sequence of piecewise constant 
approximating functions converging uniformly on X. Therefore 

J K d{Tp) =Y.p, j K d{p o i;-i) 

= '^Pi d{p oT^^) = '^pi / h„ o Ti dp. 

On the other hand, / hn d{Tp) J h d{Tp) as rt ^ cx3, while for each pair i, n the function hn o Ti 
is piecewise constant. Thus the sequence of functions {/i„ o Ti}„ converges uniformly on n to a 
function ho Ti which yields the convergence 



y^Pi J hnoTidp ^"^Pt J hoT, dp. 



The following result is a simple generalization of the well known Hutchinson Theorem [[l3[ . 



Lemma 2.2 ph{T ^,Tv) < A • pH{^^^v) for any measures G Ai. In particular contraction 
on average yields the strict ergodicity of the random map T . 



Proof. Introduce the notation 

H := {/i e C°(X, R), \h{x) - h{y)\ < p{x, y), Vx, y ^ X] 

Then 

p(T/x, Tv) = sup |y h d{Tfi) - J h d(Tv), heU 

= sup I Pi y h o Ti dfi ~ Pi y hoTi diy, /i G | . 

Consider a function h Xt ° ^'^^ each pair of points x,y ^ X we have 

< ^^K|/ior,(a;)-/ioT,(y)| 

< - JzYIp'^t, ■ p{x,y) ^ p{x,y). 



Hence h E H. Introducing another set of functions 



n-.^i^hen-. 3hen: h^ ^J2p,hoT)j , 



we can rewrite the distance between the images of the measures as foUows 

Ph{Th,Ti^) = supj^y hdp - jzj ^^'^ -hen 

Now, since C we come to the desired estimate 

Ph{Tp, Tv) < Aj7 • ph{p, i^), 
and thus the contraction on average yields the uniform contraction in the space of measures. 



Note that the exa mpl e 1.1 satisfies the conditions of Lemma 2.2, while for the example 1.2 the 
conditions of Lemma 2.2 hold only when Q < p < 1/3. Indeed, 



^p,A,^2p+hl-p) 



3 1 , 
-p + - < 1. 
2^ 2 



On the first sight it seems that the continuity of the maps Ti was not used in the proof, however 
it plays a very important role in it. One can easily construct an example when the absence of this 
property leads to the nonuniqueness of the invariant measure. 

Example 2.3 Ti{x) ~ f l[04/2](a;) + 2iil(i/2,i] (x), ^2(2;) ~ x. 



For each p G (0, 1) the random map corresponding to example 2.3 possesses exactly two ergodic 
invariant measures concentrated at points and 1 respectively. 

The contraction on average, that we assume in this section, does not prevent some of the maps 
Ti to be expanding. Therefore, despite the fact that some of the maps Ti may possess several (not 
necessary a finite number) ergodic invariant measures, the random map T under the assumptions 
of Lemma 2.2 is strictly ergodic. 



3 Expanding on average random maps 



Results obtained in the previous section are based technically on the contraction on average prop- 
erty. Now we are going to show that the opposite assumption about the expansion on average leads 
to the ideologically close result - the existence of the absolutely continuous invariant measure. 

For the sake of simplicity we shall restrict ourselves here to the analysis of piecewise C^-smooth 
maps of the unit interval [0, 1] into itself with nondegenerate expanding constants 

At. := inf |7;'(a;)| > A > 0. 

X 

It is straightforward to show that the transfer operator corresponding to the random map T in the 
space Li can be written as 

i 

where is the Perron-Frobenius operator, corresponding to the map Ti and describing the dy- 
namics of densities of measures under its action (see a detailed discussion of properties of these 
operators for example in |^). 

Denoting by Q the set of values of the index i, we consider for a given number A the sequence 
of sets 

1](")(A) ■.= {u;en: |(T>)'| > A" for a.a. x G X}, 



where 



T^x := o T^^_-^ o ...oT^^x 



is the n-th point of a realization of a trajectory of the random map T starting from the point x, 
and introduce the following regularity assumption: there exist two constants A > 1 and C < oo 
such that 

> 1 - Ce^^ (3.1) 

for each positive integer n. 

Denote by var(-) the standard one-dimensional variation o a function and by BV the space of 
functions of bounded variation equipped with the norm || • ||bv '■— var(-) -I- || • 

The following result gives the decomposition for the transfer operator under the considered 
assumptions. 



Theorem 3.1 JT^ , [J/ Let the regularity assumption {3J_) holds. Then for each pair of positive 



integers n, k the following decomposition takes place for the random map T: 
and 

var(P„,fe/i) < Const (a" var(/i) + f3''\\h\\) , (3.2) 

||g„,fc/i|| < Const\/fce-^||/i||, (3.3) 

for each function h E BV and 0<a<l</3<oo. All constants above depend on the choice of 
the map Ti, but do not depend on n and k. 

One can find in the literature dedicated to the question of the existence of absolutely continuous 
invariant measures in our setting two types of sufficient conditions for this existence. The first of 
these conditions obtained in pT[ corresponds to the strong expansion on average 

Ef<l, (3.4) 

i ^ 

while the second, described in p9| , is a weaker condition: 



One can easily show that the first of these assumption yields the second one. On the other hand, 
as we shall show there is an important difference between properties of invariant measures and 
respectively random maps under these assumptions. To explain the difference let us return to our 



regularity assumption and show that it is even more general with respect the inequality (3.5). For 
this purpose we shall need the following simple technical estimate. 

Lemma 3.1 Let be a sequence oj independent identically distributed (iid) random variables 

having exponential moments up to some positive order sq, i.e. E [e''^] < oo for all < s < sq. 
Then for each number i? > E [^i] there are constants a < 1, A < oo such that 

1 " 

V{-y^i^>R}<Aa^ 
1=1 

for each positive n. 

Proof. By the exponential Chebyshev inequality 

1 " 

> fi} < e^'^E [e'**] 

for each positive number s. For our purpose it is enough to show that the right hand side of this 
inequality decreases exponentially fast. Note that for each number x the following inequality holds 

|e"-l-a;| < el"l - 1 - |a;|. 

Indeed, this is trivial for x > (since > +x and e^^ < e^), while for x < we have < +x. 
Thus the inequality can be reduced to 

-e^ + 1 + a; < e^"" - 1 + x or e"^ + > 2, 

which is evidently correct. Therefore 

|e^«-l-sei < el^^l - 1 - |s^|. 

Assume first, that the values S,i are bounded from below. Then 



E 



1 - E[|se|] < oo 



for |s| < So, which yields the negativity of the left hand side of the previous inequality. Therefore 
there are such positive constants s £ (0, sq) and C that 

E [e^«] < 1 + E [C] + < e^EK]+Cs^ 

Thus, setting 



s = (i?-EK])/(2C), 



_ ^-(i?,-EK]^)/(4C)) 



a — e 



we get the desired estimate. Observe that s < sq by the consruction. Therefore our inequalities 
make sense only if 

R<Rq :=E[C] +2CSo- 

This means that larger values of R should be changed to Rq. To finish the proof note that if the 
random values are not bounded from below it is enough to 'cut' them from below by means of 
some constant and to apply the above argument to the result. g 



Lemma 3.2 The inequality (3^) implies the regularity assumption (3.1) 



Proof. It is enough to apply Lemma 3.1 to the sequence of iid random values := IuAt; 



On the other hand, the following result shows that the opposite statement does not hold. 



Figure 1: 

Example when the regularity condition holds while the condition (3.5) breaks down 



Example 3.3 



2x, ifO<x<^ 
ri(x) <; I - fx, tfl<x<l 



2x (mod 1) , otherwise, 
while T2{x) := Ti{x) + 1/12 (mod 1) . Graphs of these maps are shown on Fig. ^. 



Lemma 3.4 The random map T in the example 3.c for any < p < 1 satisfies the regularity 



assumption, while the condition [S.t ) breaks down. It is interesting that in this example the stronger 



assumption (3.4) holds for the second iterate T 



Proof. Both maps Ti are piecewise linear and the moduli of their derivatives take only two 
values 2 (on the intervals (0, j) and (i, 1)) and | (on the interval (i, 1)). Since both these maps 
transform the interval {j, ^) into (0, j) and on the remaining interval (i, 1) the derivatives of both 
maps is equal to 2, it follows that for any i,j G {1,2} the inequality 

|(r,r_,a;)'|>2.^ = ^>l 

holds. Thus we have checked the regularity assumption. Now observe that the derivatives of both 
maps is strictly less than 1 on the interval a; € (|, which contradicts to the condition ( |3.5| ). 

It remains to check the condition ( [3.4[ ) for the second iterate , which turns out to be a 
consequence of the fact that expanding constants for the maps (TiTjx) are not less than f > 1 
(according to the inequality above). g 



Note that since the right hand side of the inequality (|3.2| ) contains the term jS'^ with /? > 1 the 
Theorem guarantees only estimates of the type 

YaiiPnkh) < Const , 

1 — a 

which means that despite the fact that the density of the invariant measure is integrable, it might 
be not a function of bounded variation. In fact, the example 2 from Section |^ demonstrate this 
phenomenon when the parameter p belongs to (^, |). Moreover, in this example the density not 
only not a function of bounded variation, but it goes to infinity in the vicinity of the origin. 

It turns out that under a stronger assumption ( ^.4[ ) the standard Lasota-Yorke inequality is 
valid for the random map and thus the invariant density is a function of bounded variation. 

Theorem 3.2 Let for a. a. x E [0, 1] the inequality 

V^^<7<1 

^\n{x)\ -^^ 



holds. Then there are constants C, /3 < oo such that for each n G TZ^ and a function h € BV the 
Lasota- Yorke inequality holds: 

var(P^/i) < C7"var(/i) + /3||/i||, (3.6) 
from where (as usual) it follows that for each nonnegative function h £ 'L^ with \ \h\ \ = 1 the limit 

^ n— 1 

lim - P^h -. Htt, 

n— )-oo fi ^ — ^ ^ ^ 
fc=0 

exists and 

Pjvhjv = /ijT, var(ft,jT) < Const , 

while which respect to the absolutely continuous invariant measure fj-j; with the density hjT the 
correlations decay exponentially. 

Now we are able to finish the analysis of the example 2 from Section Observe that for p > \ 
we have 



which yields the condition ( |3.5| ) and, hence, our regularity assumption. Therefore for p > 1/2 
there exists an absolutely continuous T-invariant measure. On the other hand, for | < p < 1 the 
condition (3.4): 



holds for a. a. x £ [0, 1]. Thus we can apply Theorem 3.2, whereis the boundedness of the density 
of the invariant measures follows. 



4 Perron-Frobenius spectrum (PF-spectrum) 

In the previous sections we have restricted the analysis of statistical features of random maps to the 
properties of their invariant measures and more specifically Sinai-Bowen-Ruelle (SBR measures). 
From a more general point of view the SBR measure is the eigenfunction of the Perron-Frobenius 
(transfer) operator of our random map in a suitable Banach space corresponding to the leading 
eigenvalue (1). Therefore it is very natural to extend the analysis of the dynamics to the complete 
spectrum of this operator. 

It is worth note that the interest to the PF-spectrum is based to a large extent on the fact 
that the subleading elements of the spectrum define the rate of mixing (convergence to the SBR 
measure, correlation decay, etc.). 

4.1 Definition of the PF-spectrum 

Let us start with the short description of objects related to the notion of the spectrum which 
we shall need further. Let P : S — *■ S be a bounded linear operator in a complex Banach space 
(B, IHIg). As usual we denote by spg(P) its spectrum, which is defined as a complement to the 
set of regular elements, i.e. to the points z € C such that the resolvent {zl — P)~^ of the operator 
is defined in the entire space and hence is bounded. The maximal (on modulus) element of the 
spectrum is called the spectral radius: 

Pb{P) :=sup{|z| : z G spe(P)}. 
As it is well known ^ the spectral radius can be calculated by the following formula: 

Pb{P)^ lim ||P"||^/". 

n — >QO 

Browder ||^ introduced the notion of essential spectrum cssSPb(P) of a bounded linear operator 
P as a union of the elements of the spectrum z E spg(P) such that at least one of the following 
properties holds: 



1. The region of values of the operator z/ — P is not bounded in B. 



2- U„>oker((z/ — P)") is infinite dimensionaL 

3. 2 is a hmit point of the spectrum spg(P). 

Outside of the essential spectrum only a countable number of isolated eigenvalues of the operator 
P may occur. Naturally the essential spectral radius jOB^ess(P) of the operator P is defined as the 
minimal nonnegative number such that all elements of the spectrum spg(P) outside of the disc 
{z e C : \z\ < pe,oss(P)} are isolated eigenvalues of finite multiplicity. It turns out that the 
essential spectral radius can be calculated by a formula similar to the one for the usual spectral 
radius: 

PB.css(P)-,lim ||P"||^';,",,, (4.1) 

but for the specialized seminorm: 

||P||g := inf{||P — : K : B ^ B - a compact operator}. 
Clearly for each e > the set 

sp(P)n{ze C : |z| >pe,css(P)+e} 

consists of a finite number eigenvalues Ai, A2, . . . , \M(e) of the operator P. Schematically on the 
complex plane the spectrum can be represented as a disk of radius pe,css(P) centered at the origin 
(describing the essential spectrum) and (not more than countable) collection of points between 
this disk and the circle of radius pb(P) also centered at the origin. 



4.2 Spectrum for the case of contracting on average random maps 

In Section |^ it was shown that a contracting on average random map T possesses the only one 
invariant measure /i^r to which the sequence of iterations {T /x} converges (exponentially fast in the 
Hutchinson metrics) for each probabilistic initial measure jjl. From the point of view of dynamical 
system theory the following question after this is the analysis of the spectrum of this convergence. 
The problem here is that typically the limit measure /x-y is not absolutely continuous, which rules 
out the description of the transfer operator P^r in a space of a reasonably 'good' functions, for 
example, in the space of functions of bounded variation. To overcome this difficulty we shall study 
the action of the operator P^r in a much larger space of generalized functions equipped with a 
norm induced by the Hutchinson metrics. 

Let (X, p) be a d-dimensional smooth manifold with a finite collection of continuous maps {T^} 
having bounded Lipschitz constants At^ from X into itself, and a collection of probabilities {pi], 
defining the random map T. Remind that the contraction on average means that 

Ajt:=^P,At, < 1. 

i 

Before to define our space of generalized functions we need first to define the class of test- 
functions if : X ^ M^. For this purpose we introduce the following functionals: 

rr I \ |y?(x) - ^){y)\ 

^"('^) - sup , ,,^ ' 

x.y^X, p{x,y)<v H x^^tUI 

where \'f\oo '■= esssup|(p|, and the constant v £ (0, 1]. The first of these functionals is the Holder 
constant with the exponent a < 1, while the second one for each finite nonnegative value of the 
parameter a is the norm in the Banach space of a-H61der functions on X, which we shall denote 
by C". Without the loss of generality we shall assume that the diameter of the phase space 
sup^ y^x Pi^^y) ^ 1- Note that the another restriction p{x,y) < v is introduced only to be able 
to work with the exponential map on general smooth manifolds: in the case of a flat torus this 
restriction can be omitted (or simply one can set u = \). 



Consider now the space of generalized functions on X with the norm defined in terms of the 
test-functions from the space C" : 

\M,. ■■= sup / h<p. 

V„(v)<li 

The proof that the functional \ \-\\^^^ is indeed a norm in this space is standard and we leave it for 
the reader. 

Denote by J^ce the closure of the set of bounded in the norm 1 1 • 1 1 generalized functions from 

jr. 

Lemma 4.1 J^p C To. for any numbers < a < /3 < 1 and each function G C the inequality 
V/si^) < Vai^p) holds. 



Proof. Indeed, 



Hx)-(fi{y)\ \ip{x) - ip{y)\ 



Extending the standard definition of the Pcrron-Probeniusa operator to the action in the space 
of generalized functions we get the representation 



j Pj^h-ip = j h-Y^PiipoTi=: j h-{ipoT) 



for each test-function from ^p e C". 

Let us fix some constants 0<a</3<l, 0<a<oo, whose exact values we shall define later. 
Our first aim is to derive a version of the Lasota-Yorke inequality for the action in the space Ta- 
For g e [0, 1] define a function 



which wc shall need in this derivation. 

Introduce additionally the notation G = {Gi,Pi} for the random map defined by the collection 
of maps Gi and the distribution {pi}. 

Lemma 4.2 The superposition of any pair of random maps G = {Gi,pi\ and G = {G^,p'j} 
acting on the same manifold {X, p) and satisfying the Lipschitz condition is again the random map 
g' oG:= {G'i o Gj,p'iPj} and A^'^^(g) < A^/(g) • A-^{q) for each q e [0, 1]. 

Proof. Indeed, 



Theorem 4.1 For each number k> 2 and for any h € J^a and n € !S+ the Lasota- Yorke inequality 
holds: 

< KA^{a)\\h\\^^^ + Const • - 2)-V« \\h\\^^^ . (4.2) 



Proof. We start from the proof of the foUowing two inequahties: 



II I 11(3) - " ••(p) 

\\Pt^\\^ ^ < KAj.{a)\\h\\^^^ + Const • (k- 2)"^/" 
By the definition of the Holder constant we have: 

\(p{Tix) - (p{Tiy)\ p'^{TiX,T,y) \(p{Tix) - (p{Tiy)\ 



(4.3) 
(4.4) 



p"{x,y) p"{T,x,T,y) 



Hence, 
Thus, 

On the other hand, since 
then 

Therefore for each (3 £ (0, 1] we have 



p°'{x,y) ^ 



oo S \'fi\< 



which imphes the incquahty (4.3) for each /? £ (0, 1] and a function h e Tp. 
To prove the second inequahty ( |4.4| ) we need more dehcate estimates. 
Introduce the following notation: 

Bs{x) -.^{yeX: p{x, y) < 5} Bg := {C G %X : |^| < S}, 

i.e. Bs{x) is the ball of radius 6 centered at the point x in the space X, while Bs is the ball 
of radius 6 centered at the origin in the tangent space TxX. For each point x G X consider the 
exponential map 

:= exp^ : Bg C %X X. 

Choosing 5 > small enough we always can assume that 'i'xBs C By{x). 

Denoting by m the Lebesgue measure on X, we introduce the following smoothing operator 



}sv{x) := 



m{^.xBs) 



where J^x is the Jacobian of the map '^x- 

Let us estimate the Holder constant of the function Qgifi'- 



IQsifix) ~ Qsipiy)\ < 



J^x{z)dz- I J^!y{z)dz 

Bs JBs 



L </5(*yZ) • J4'y(2) dz 



1 



/ J-^x{z)dz 



IBs JB. 

^ vT'-i'T' n-l I liT,-i 

<2 supJ*^ • |(p|oo 



/s, J-^xiz)dz- J^^ J'ify{z)dz 

(fii^xz) ■ J-^x{z) dz- fi'i'yz) ■ J'i'yiz) dz 
JBs 

\Bs\^-'^yBs\ + \^-'^yBs\Bs\ 



< ^Ci\<f\oop{x,y), 



where the constant Ci depends only on the properties of the manifold X. 

Now we estimate how much the operator Qg differs from the identical operator: 



\Qsip{x) - (p{x)\ < 



where the constant C2 also depends only on the properties of the manifold X. 
On the other hand, 



Thus, 



< 2A^(a)i/„(^)+C72^"i?a(¥')- 



Gathering the obtained estimates and using Lemma 4.1 



we come to 





< sup / 







sup Vc,ly^p,{ipoT,-QsiipoT,)))] 
.M<i \i J 

\\h\\^^^+{2A^{a)+C,6'^)\\h\\^^^. 



Va 

Ci 



The refo re choosing the value of the parameter S such small that 6*25" = k — 2, we get the inequal- 
ity 0. 



Now according to Lemma 4.2 and above inequalities we get 

= ||P5T"/i||^ ^ < KA^..(a) + ConstJ-i 



< AC AMI I 1 1 + Const • (k - 2)"!/" 



which finishes the proof of the inequality (4.2) 



Lemma 4.3 The function Ajt(-) is convex, takes values strictly less than 1 in the interval (0, 1], 
and under the condition 



Ha 



Pi At 



< 1 



(4.5) 



its unique point of minima either lies inside of this interval, or is larger than 1 otherwise. 



Proof. By the definition of the contraction on average we have Ajt(1) = Ajt < 1. On the other 
hand, A^(0) — 1, and 




A typical behavior of the function A7p{q). 



Figure 2: 



Thus the function A^(<7) is strictly convex and for any q G (0, 1] it is less than 1. Let us show that 
this function strictly decreases at 0. Indeed, 

i 

On the other hand, from the contraction on average 

^P,At, < 1 

i 

using the convexity of the logarithmic function, we get: 

^PilnATi < Inl = 0, 



which proves that ^Ajt(O) < 0. It remains to check the last statement. 

^A5.(l) = ^^At, • In At, ^ ^lnA?,f-- = In (l] ^^rf 



Thus due to the inequality (4.5) the unique (due to the strict convexity of the function Ajt(-)) 
solution of the equation j-A^(q) — (the point of minima of the function A^(-)) belongs to the 



interval (0, 1). On the other hand, if the inequality (4.5) does not hold this solution is greater than 
1. ■ 



A typical behavior of the function ATp{q) is shown Fig. ^. Denote by q the value of the parameter 
q corresponding to the unique (due to the strict convexity) minima of the function A^(q). The 
value q is positive since ^Ajt(O) < 0. The position of q with respect to 1 is defined by the sign of 
the derivative of the function Aj^iq) at 1 which might be both positive and negative. For example, 
if all the maps are contractive then this sign is negative and the function Aj^^q) strictly decreases 
on the interval [0, 1]. On the other hand, the map in the example 1.2 satisfies the condition ( |4.5| ) if 
< p < 1/5 and in this case q = ^ log2(l/_p — 1). We describe the properties of the function A^(g) 
in such detail because the fact that it can grow in the vicinity of the point 1 plays an important 
role in further calculations. 

Lemma 4.4 The unit disk in the strong norm \\-\\^ ^ is a compact set in the weak norm \\'\\^^^ 

The proof of this statement follows immediately from standard results on the enclosure of the 
spaces of Holder functions. 



Above statements imply by the lonescu-Tulchea and Marinescu Theorem |14 the quasicom- 
pactness of the operator P^r, and the vahdity of the following based on Nussbauni Theorem pO| 
estimate of its essential spectral radius. 

Lemma 4.5 The essential special radius of the operator P^r : J-a belongs to the disk 

or radius ATp{a) centered at zero. 

Note that the estimates leading to the Lasota-Yorke type inequalities depend sensitively on the 
choice of the value of the parameter a. Thus it is reasonable to choose the value of a which yields the 
smallest (and hence the best) available estimate of the essential spectral radius. Normally (compare 
to W) this value is equal to 1 which is unavailable since we consider only Holder continuous test 



functions. However in our case Lemma L3 shows that under the condition ( |4.5| ) the optimal value 
of a may be strictly less than 1. 



An immediate corollary to Lemma L5 is the existence of a constant 7 g [AjT(a), 1) such 
that the set sp(Pj^) \ {\z\ < 7} consists of a finite number of periferal eigenvalues ri,...,rjv 
of finite multiplicity. Denote by Pi, . . . , Pn the corresponding spectral projectors and set P := 
^ Si^i ^3- Then the rank(Pj) < 00, P^Pj — rjPj {j — 1, . . . , N), and the spectral radius of 
the operator P^P does not exceed 7. 

Besides, 

• If |r| = 1 the operator 

1 ^ "'^ ^ \ n — l ^ ( P f 

r / 1 0, otherwise 

k=o j=i fc=o ^ 

is well defined in the ||-||^ ^-norm. In particular, P, = P^^, and since / |Prj/| < / I/I for all 
/ e C^{X,1R.^), then operators Pj can be extended continuously to the entire space L^. 

• For any function / G PjJ'a there is a finite Borel signed measure /i/ on AT such that (/, (p) = 
J (pdfif for all ifi e C^(A, M^). ri := 1 e sp(P^), fj, := ^ip^i - a positive measure, /i(A) = 
m(X), and all signed measures /i/ are absolutely continuous with respect to /i. 

One can interpret these statements as follows: for f,ip€ C^(X, R^) and |r| = 1, 

-J n — 1 -. n— 1 ^ 

(P./, ^) = hm - J2 r-'iP^f, ^) - lim - ^ ^ o p' • /. (4.7) 

fc=0 fc=0 

Hence, {Prf,ip) < |<p|oo • / |/|, so Prf can be extended by continuity to a continuous linear 
functional on C"(A, H^), and by Riss Theorem there is a measure fJ-p^f, such that {Prf, f) = 
J (pdfip^f. If r = 1 and f,(p>0 then J (pd^p-^f = (Pi/, ip) > and (Pi/, 1) = / / according 
to ( [4.7[ ). /ipii is a positive measure and ?' = 1 is an eigenvalue of the operator P^r. Finally, 



it follows from (4.7) that for any 1^ > we have 



|(P,/,^)| < limi|] /^op'-l/l < I/loo (Pi 1, V') = I/loo / ^■ 
Besides, /ipj/ is absolutely continuous with respect to fj.. 

It remains to show that PjJ-a = Vj := Pj(C^(A, H^)). Since Vj C PjJF„ and they are 
finite dimensional linear subspaces in J-'a, then this statement immediately follows from the 
denseness of the space 

• T = 11, since Jipd{T 11) — J ip o T dfi = limn^oo ^ X]fe=o / 'z' ° ^ ^ — Jf^lJ. for all 

ip e ci(x,iRi). 

• If r = 1 is a simple eigenvalue and there are no other eigenvalues equal to 1 on modulus, then 
Pi/ = (/, 1) -Pil for all f eTa- This follows immediately from the fact that (Pi/, 1) = (/, 1) 



(see (^ 



c 



1/2 



1 - c 



Figure 3: 

Example of a random map with a nontrivial isolated eigenvalue 
• is a SBR measure, since for each </? €E C-'^(X, H"'^) 



lim / ipdi — r m ) = lim — / (p o t ' — {Pi l,Lp) — [ 

n^oo I \ n ^ — ' / ri^oo n ^ — ' / / 



(fdfi 



k=0 ' k=0 • 

So far the only isolated eigenvalue that we were able to identify explicitly was the unit eigenvalue 



(the leading one). Let us discuss the following argument. According to Lemma 4.5 the essen tial 
spectral radius of the operator Pjt is not larger than A^ia)- On the other hand, the Lemma 2.2 
garantees the convergence to the limit measure with the rate at least 1 > > AijT{a) (by 
Lemma 4.3). Thus if there is a function / e such that Pjt/ = A^r/, then A^r is an isolated 



eigenvalue. 

Example 4.6 Ti{x) := l/2 + sign(x- 1/2) • [1/2- |2sign(x- 1/2) • (a;- 1/2) - 1/2|], where sign(x) 
is the sign of the number x, and T2{x) :— x/2. 

The maps Ti are shown on Fig. ^. Note that this random map is contractive on average when 

< p < 1/3. By the construction the map Ti has three unstable fixed points c — 1/6, 1/2 and 

1 — c = 5/6 and the point 1/2 is the only fixed point of the contractive map T2. Besides, the 
trajectories of the random map starting at points c and 1 — c are symmetric and consist of a 
countable number of points {2~^ — 2~"c}„ and {2^^ + 2^"(1 — c)}„ respectively. Consider a family 
of generalized functions 

n n 

fa{x) := ^flfc • l2-i_2-fcc(a;) - ^ak- l2-i+2-''(i-c)(a;), 

fe=0 fc=0 

parametrized by the sequence of coefficients d — {ak\k- Due to the previous remark this family 
is invariant with respect to the action of P^r and our aim is to find such a sequence d that 



Lemma 4.7 In the example \4 (\ for each < p < 1/5 there is a sequence a — a{p) such that 
PjT/a=A^/a. However WfaW^^^ ^ 00 for any a G (5 log2(l/_p - 1), 1). 

Proof. Observe that A^r = (1 + 3p)/2 < 1. Therefore rewriting the eigenvalue relation in terms 
of the weights ak we obtain the following recurrent relations: 

((1 + 3p)/2) • flo =pao+pai 

((1 + 3p)/2) • flfc ^ {1 - p)ak-i + pttk+i 

for each fc > 1. Solving these equations with respect to the variables with higher indices we get: 

1+p l + 3p 1-p 

ai = — — flo, Qk+i = — flfc ak-i lor fc > 1. 

2v 2v V 



/ l+3p 1-p \ 

The eigenvalues of the matrix A ~ i "^p (f j ' controlhng the growth of the coefficients, are 

equal to 2 and (1 —p)/{2p) respectively. Now since for < p < f/5 the second eigenvalue is larger 
and since ai/ao — (1 +p)/{2p) > 2 we deduce that the constants grow as ((f — p)/{2p))^ as 
k — > oo. 

It remains to estimate the ||-||^ ^-norm of the generalized function fa'. 

oo 
k=0 

if and only if ((1 - p)/(2p))2-'=" < 1. On the other hand, (1 - p)/{2p) > 2 for < p < 1/5 which 
contradicts to the convergence. | 



4.3 Spectrum for the case of expanding on average random maps 

Let Ti for each i be a map from the unit interval into itself and let the following condition holds: 

where the supremum is taken over all points x Cz X where the derivatives of maps T.^ are well 
defined. Then, according to Theorem |3.2|, the Lasota-Yorke inequality is valid: 



var(P^/i) < C7^var(/i) + /3||/i||. 

Therefore all known results on the spectral properties of the Perron-Frobenius operator obtained 
for piecewise expanding maps based on a similar inequalities remain valid as well (see for example 
the detailed discussion in 



As we already mentioned the condition (3.5) does not imply the inequality of Lasota-Yorke type 
and moreover there are examples when under the condition (|3.5| ) there is no exponential correlation 
decay. Therefore the question how to extend the description of the spectrum to this case remains 
open. 

4.4 Stochastic stabihty 

In this section we shall study random perturbations of random maps under consideration. Since 



under the condition (3.4) perturbations of expanding on average random maps can be considered 
exactly as in the case of deterministic piecewise expanding maps (see ^, ^ and general discussion 
in ||l|) we shall restrict ourselves only to the case of contracting on average systems. 

As usual under the randomly perturbed system we shall mean the superposition of the original 
system and a Markov process acting on the same phase space and defined by the family of transition 
operators Qe (here e stands for the 'size' of perturbation). 

To simplify the calculations we shall start from the case X = Tor'' and then shall explain how 
the corresponding arguments should be changed in the case of a general smooth manifold. 

Consider two families of operators: integral operators Qe : JF„ ^ and the dual ones Q* : 
Ci^Ci: 

Qef{x) := y qe{z,x)f{z)dz, 

Q*eV{x) := y" qe{x,z)tp{z)dz 

with the family of nonnegative kernels qd', ■), with respect to which we shall assume that for some 
1 < Af < oo the following conditions hold: 

qeix,y)dy = l, qe{x,y)^0 yp{x,y) > e, (4.8) 



\qeix,z) '-qdy,z + y~x)\ dz < Mp{x,y), (4.9) 
We start the analysis of the operator Ql from the following simple estimates. 



Lemma 4.8 For each e C" we have 

\Q*e<pU < IV'loo, (4.10) 

\Ql^-v\oo<£'^H^{^)- (4.11) 
Proof. The proof is straightforward: 

\Q*e^{x)\ =1 j qe{x,z)ip{z)dz\< j qe{x,z)\ip{z)\dz <\ip\oo. 

\Ql'f{x) - ip{x)\ =1 j qe{x,z)^p{z)dz- ip{x)\ 
< j qs{x,z)\ip{z) - ip{x)\dz 

since only the points z £ B^ix) should be taken into account. g 

Now let us estimate the norm of the operator Q^. 
Lemma 4.9 Let Mi(e) := max{2e«/2, Me^i-")/^}. Then \\Qe\\^^^ < l + Mi{e). 

Proof. Our aim is to show that is a valid test function and to estimate the values of VaiQ^i^)- 
We already estimated the supremum norm of the function To get the estimate of the Holder 
constant we consider two different situations: when the points x, y are close, i.e. p{x, y) < \fe, and 
the opposite case when they are far apart. In the first case we proceed as follows: 

|Q*<p(a;) - Ql^{y)\ = \ J qe{x, z)ip{z) dz - J qe{y, z)ip{z) dz\ 

= I j (le{x, z)^p{z) dz- j qe{y, z + y- x)ip{z + y-x) dz\ 

< I j qe{x,z)ip{z)dz - J qs{x,z)ip{z + y-x)dz\ 

+ I j Qeix, z)ip{z + y-x)dz- j qe{y, z + y- x)if{z + y - x) dz\ 

< j qs{x,z)\ip{z) - (p{z + y - x)\dz 

+ / \qe{x,z) - q^{y,z + y - x)\ ■ \ip{z + y - x)\dz 



< p"{x, y)H„{^) + Mp{x, yM^ 

< H^{ip) + M£(i-")/Vu] p"(a;, y). 

In the opposite case, when p{x, y) > y/e we shall proceed in a different way: 

\Q*Mx) - Q*My)\ < \Q*Mx) - v{x)\ + \Q*My) - >fi{y)\ + Mx) - <p{y)\ 
<2\Q*<f-ip\^+p"{x,y)H„{<f) 

<2£"Jf„(^)+p"(x,y)F„(^) 

< (i + 2£"/V(a;,y)ir«(<^). 



Hence we have 



< (1 + 2e°'/'^)Ha,{>f) + M£(i-")/2 

< (1 + max{2£"/2, M£(i-")/2})y^(<^). 



Thus, setting Mi(e) max{2e"/2, Me(i-")/2}, we get 

IIQe/IL , = sup Qef-f= sup f-Ql(p 

< sup yalQ*"^)- sup 



.(V)<1 K.(V')5 

< (1 + Mi(e)) 



Observe that in several places we used the estimates of the supremum norms obtained in Lemma 4.8. 



Lemma 4.10 Let G : J-q Ta he a linear operator, and let G* : — > be dual to it, i.e. 
jGf -ip^ J f -G*^. Then for all f e 

< ( sup 1^„(GV)| .11/11 

Proof. Indeed, 

\\Gf\\ = sup [g/-^^ sup //-GV 

Vfi{f)<lJ Vff{v)<lJ 

< I sup T4(GV) ) • sup //-v. 



Lemma 4.11 



\\\Qe 



WQe 



sup ||Qe/-/|L 



as e — > 0. 



Proof. Applying Lemma 4.10 to the operator G = Qe — 1 we get that the sufficient condition of 
the validity of the desired statement is the convergence of 

sup Va{Ql(p - (p) ^ 

Vi,{v)<l 

as e — > 0. Let us prove this convergence. Observe that since fi > a and Lp e we can get a 



stronger estimate compare to Lemma 4.8 



\Q*eV{x) - 'P{x)\ =1 j qe{x,z)ip{z)dz- Lp{x)\ 
< / qe{x, z)\lp{z) - Lp{x)\ dz 



Applying now estimates similar to ones used in the proof of Lemma |4.9| and taking into account 
that we consider more smooth test-functions ip £ in the case p{x, y) < e we get 

\{Q*M^) - ^i^)) - (QMy) - ^iy))\ 

< \Q*Mx)~Q*My)\ + Mx)^p{y)\ 

< /(x, y)Hp{p=) + Mp{x, + P^{x, y)H[3{p>) 

< [2eP-"Hp{p)+Me'-"\cp\oo] p"(a;,y). 



While in the opposite case, when y) > e, using the same argument as in the proof of Lemma 4.9 
we get 



<2e^-^Hp{^)p'-{x,y). 

Hence for Lp S 

which yields the following estimate 

< Se'^-'^H^iifi) + M£i-"|v?U < (3 + Me^-''^)e^-"Va,{^) ^ 

as e ^ 0. 



The properties of the transition operator obtained above together with Theorem 4.1 under 
the additional assumption that Ajt(q;) < 1/2 make it possible to use results about the spectral 
stability of transfer operators satisfying Lasota-Yorke type inequalities ||l^ and to obtain the 
following stability result. 

Theorem 4.2 Let the conditions U-^ satisfied and let A-^ia) < 1/2 for some a € (0, 1). 
Then all elements of the spectrum spjp^(Pjt) outside of the disk of radius Ajria) are stochastically 
stable and the corresponding eigenprojectors of the perturbed system converge to the genuine ones. 



Proof. First let us show that the transfer operator for the stochastically perturbed system 
satisfies a Lasota-Yorke type inequality. A straightforward calculation shows that this operator is 
equal to QeP^. Combining the results of Lemma 4.9 and Theorem 4.1 we get for any h € that 



< il + Mi{e))\\Pj;h 



< {1 + Miie))KAj;{a) 



Const • (k - 2) 



-I/O 



'(/3) 



Therefore, if A^(a) < 2 for some < a < 1, then the number 7 := (1 + A/i(e))KA^(a) < 1 for the 



value of K > 2 guaranteed by Theorem [O 
stability result in 



Since all other assumptions of the abstract spectral 
|15| were alredy checked during the analysis of our Banach spaces of generalized 
functions, we come to the desired statement. g 



To proceed further we need to generalize the notion of the periodic turning point, well known 
in the one-dimensional dynamics. Namely, a point x Cz X is called the periodic turning point for 
the map T : A" — s- X if T"x — x for some n G and the derivative of the map T is not well 
defined at the point x. 

Definition 4.1 A point x E X is called the periodic turning point for the random map T : X ~* X 
if there is a finite collection of indices ii,i2, ■ ■ ■ ,ik such that the point x is the periodic turning 
point for the deterministic map ^ o Ti^ o . . . T,;^ . 

For example, the point x — 1/2 is the periodic turning point for the random map in the example 
2 for any nontrivial distribution (0 < p < 1). 

Theorem 4.3 The assumption A^(a) < 1/2 can be replaced by the following: either all the maps 
Ti are bijective and -differentiable, or they are piecewise -differentiable and have no periodic 
turning points. Then all isolated eigenvalues are stochastically stable. 



Proof. The key idea here is to consider another representation of the perturbed operator 
(QeP^)" = (5g(„)P^ and to show that the new operator satisfies the same assumptions 

( |4.8| , |49| ) as the operator Qe (except that the value of £{n) might be in 1/ min{ATi} times larger). 
Therefore choosing n large enough we can always get A^(a) < 1/2. This idea was first applied in 



the case of piecewise expanding maps in g, g| and then in the case of hyperbolic maps in Q . If 
all maps T; are bijective this can be done by a simple change of variables, while in the second case 
the construction is more involved but is completely similar to the one in |S, p[ . g 



Let us show now what should be changed in the case of a general smooth manifold. Since 
locally in a neighborhood of a point x & X one can introduce local coordinates by means of the 
exponential map '^x^ the tangent linear space T^X can be isometrically mapped into M''. Let 
> be a number such that for each point x & X the ball (v the metrics p) of radius v centered 
at this point belongs to the domain of values of the exponential map ■ Note that we already 
have introduced the restriction on the distance between the points in the definition of the Holder 
constant needing to be content with the domain of definition of the exponential map. In fact, the 
first difference appears only in the analysis of random perturbations, in particular, the condition 
( |4.9| ) should be rewritten as 

j \q,{x, y) - q,{-^x{-^-\x) + t),-^y{-^-\y) + t))\ dy < p{x, ^.(^-^(x) + t))M, 

where t eR^ and \t\ < v. 

Assuming now that e < v and replacing the expressions of type z + y — x to 

{^-\z) + ^-\y)~^-\x)), 

we obtain the same estimates as in the flat case (when X is the unit torus). Therefore all results 
of this section remain valid for the case of a general smooth manifold. 



4.5 Finite rank approximations 

Let us discuss now finite dimensional approximations of transfer operators. Again due to the same 
reason as in the previous section we shall restrict the analysis to the case of contracting on average 
random maps. 

Let {/S.i]i be a finite partitions of the phase space X into domains (cells) of diameter not 
larger than (5 > 0. For a point x £ X hy A^, we denote the element of the partition containing it. 
Under these notation the so called Ulam approximation can be described as an operator 



Qsf{x) 




Note that this operator is selfdual, i.e. Qs = Q%. One can easily check also that the dimension of 
the space Qs^a coincides with the number of elements in the Ulam partition. 

Lemma 4.12 WQsW = oo. 



Proof. Let a point yo G X belongs to the boundary between two elements of the Ulam partition, 
and let the points y{e) and y'{e) belong to neighboring elements of the partition both on the 
distance e from yo- For the function 

f,{x) := ly{e)ix) + ly'(e)ix), 

where Ij, means the 5-function at the point y, the following inequalities hold: 

IIAII^^^ < Const e", 
IIQ^MI,^, > Const >0. 

The first of these inequalities follows from the definition of the norm 1 1 • 1 1 ^ ^ , while the second 
one is a consequence of the fact that the function Qsfe is the characteristic function of the 
union of two neighboring elements of the partition containing the points ?/(e) and y'{e). Thus, 







5/12 1/2 



1 



Figure 4: 

Counterexample for the original Ulam construction. 

This result shows that the original Ulam approximation scheme cannot be immediately applied 
to the spectral analysis of random maps. 

What is still possible is that the leading eigenfunction - SBR measure may be stable for the 
class of maps we consider (the above example does not contradict to this - the SBR measure 
is preserved). I believe there should be very deep reasons explaining the stability of the leading 
eigenfunction, while all others are not stable, however presently we do not have the adequate 
explanation. In the literature (see, for example, and further references therein) the stability 
of the SBR measure is proven for the class of piecewise expanding maps. Moreover numerous 
numerical studies confirm this stability for a much broader class of dynamical systems. To the best 
of our knowledge the following simple example of a one-dimensional discontinuous map represent 
the first counterexample to the original Ulam hypothesis. 

Lemma 4.13 The map 

ff + i ^fO<x<^ 
Tx := < -2x + I if <x <} 
I f ^ i otherwise. 

from the unit interval into itself is uniquely ergodic, but the leading eigenvector of the Ulam ap- 
proximation Hi/jjFt does not converge weakly to the only T -invariant measure. 

Observe that the situation when the SBR measure becomes unstable with respect to Ulam 
scheme is indeed very exotic and the map in the corresponding example is not only discontinuous 
(see Fig. ^ , but this discontinuity occurs in a periodic turning point (compare to instability results 
about general random perturbations in [||). 

Proof. Denote by u'-"-' the normalized leading eigenvector of the matrix Hi/^Pt- A straightfor- 
ward calculation shows that for each n € 7Z^ all entries of the vector are zeros except the 
first entry, which is equal to 1/3 and the (n -I- l)-th one, which is equal to 2/3. Compare this to 
the only invariant measure of the map T ~ the unit mass at the point 1/2. g 

To overcome this difficulty we consider another 'smoothed' approximation scheme. In each 
element of the partition {A^} (of diametr < S) we fix an arbitrary point (its 'center') Xi e A^. 
Now for a given smooth enough kernel qd', •) satisfying the assumptions from the previous section 
we define the following finite dimensional operator: 

Qe,sf{x) ■.^^l^^{x) j qe{z,Xi)f{z) dz. 
Observe that the dual operator is equal to: 



Q*e.5'^{x) — ^qs{x,Xi) / ip{z) dz. 



Indeed, 

Qe,sf{x) ■ (p{x) dx = ^Ai{x) / qe{z,Xi)f{z) dz ■ ip{x) dx 



= j f{z) ^ qe{z, Xi) ip{x) dx dz = J f{z) ■ Ql^sip{z) dz. 



Lemma 4.14 Let additionally to the assumptions ( (.8, ^.i ) for any points x,y,z€X the inequal- 
ity 

\qe{x,y) - qeix,z)\ + \q,{y , x) ~ q,{z , x)\ < Me-'^-^ p{y, z) . (4.12) 

holds. Then 

WQef-Qe.sfW, , < SMe-'^-ljl-" 



Proof. 

J {Qef - Qe,sf) • = J f- {Q> - QUV)- 

Denote 

^{x) ■■= Q*eV{x) - Q*e,SVix) = ^2 ^) ~ 1^^^^ X,))ip{z) dz. 

, J A. 

Then 

l^'loo < I'/'loo • sup^ / \qe(x,z) - q^{xi,x)\ dz 
^ , Ja. 

<|(p|oo-sup^ / {\qs{x,z) - qs{x,Xi)\ + \qs{x,Xi) ~ qe{xi,x)\) dz 

^ j J Ai 



Ai 



since x,z G Ai and, hence, Tiiax{p[z, Xi), p{x, Xi)} < S. 
Let us estimate Ha{^). If p{x,y) < S then 

mx) - < \Q*Mx) - Q*My)\ + \QtM^) - Q:My)\ 

< J \qe{x,z) - q^{y,z)\ ■ \if{z)\ dz + Y^\qe{x,Xi) - qe{y,x,)\ 

< Me~''~^p{x,yMoo + Me-'''^pix,y) ■ \^\oo - 2Me-''-^d^"'' p"{x,y) ■ \ipU 
Otherwise if p{x, y) > S we apply another estimate 

mx) - $(y)| < 2|$U < 4M£-'^-i5|(p|oo < 4M£-'^-V"(a:,y) • <5'"" • kU- 

Thus, 
and hence 

which yields the desired statement. 



Theorem 4.4 Let the family of kernels {q^{-,-)} satisfies the conditions ( |^.4 , Then 

\\\Qe,s - 1||| < (5M + 3 + Me'^-^)ie-'^-^5^-" + e^-") -> as e''^~^5^~" + e^-" 0. 

Hence for the case A-jT(a) < 1/2 the isolated eigenvalues and the corresponding eigenprojectors of 
the operator are stable with respect to the considered approximation. 



Proof. According to Lemmas |4.9| and [4.14 



WQe.sW^^, < IIQell,^, + \\Qe,5 " < 1 + Ml{e) + hMe 



-d— 1 — a 



which proves the first statement. 



Similarly but using Lemma 4.11 instead of Lemma [4.9| , we get 

< (5M£-'^-i(5i-" + (3 + Me^-^)e^^-°') ■ Va{(p), 
which finishes the proof. g 

In fact the finite dimensional approximation defined by the two-parameter family of operators 
{QE,i5}e,i5 one can consider as a smoothed version of the original Ulam construction, which corre- 
sponds to the case e = 0. Observe that in our approximations the relation between the parameters 
is completely different - it is necessary that e ^ S. 

We consider also another (seeming more natural) finite rank approximation scheme. Denote 
by the pure Ulam approximation operator corresponding to the partition into domains {A^} 
whose diameters do not exceed S: 

where A^^ stands for the element of the partition containing the point x. Note that this operator is 
self adjoint. We shall approximate our transfer operator P^r by H^Q^P^. To study the properties 
of this approximation we need as usual to analyze properties of the adjoint operator, i.e. of the 
operator 



Lemma 4.15 \\Q, - IlsQ,\\^^^ < {3 + 2M)e-'^-i(5"(i-°) as e^d-igad-a) _^ q 
Proof. Denote 

<^{x) := Q*ip{x) - QtU*sip{x) = J Qeix, z) {^{z) - |^ </j(s) ds^ dz 

= j qeix, ^) 1^ i'Pi^) - '/'(^)) dz. 

Since (p € C" and the diameter of the elements of the partition does not exceed S, we have 

|$|oo < <5"i/aMsup fq,{x,z) dz^S'^H^i^). 

X J 

Now we are going to estimate the Holder constant of the function $. which we shall do in two 
steps. First, we consider the case when p{x,y) < (5": 

|$(x) - <i>(y)| < \Q*Mx) - Q*My)\ + \Q:n*M^) - Q:u*My)\ 



< J \qe{x,z)-qe{y,z)\-\ip\oo dz + J \qsix,z)-qeiy,z)\- \ip{s)\ ds dz 

< 2Me-''-^p{x,y)\^\^ < 2Me~^~^S"^'-^'>\p\ooP''{x,y). 
In the opposite case, when r{x, y) > we use a different estimate: 

< 2|<i>U < 2ril,(^) < 25"^^'"'>H^{p)p''{x,y). 

Thus 



as e-rf-ij^Ci-") 0. 



Observe that the rate of convergence in this approximation is lower compare to the previous 
one, however the numerical application of the 2ncl scheme is more straitforward. 



Corollary 4.16 Again as in Theorem ^.5 and due to the same reason the assumption A-pp^a) < 1/2 



can be replaced by either the bijectivity of the maps Ti or the absence of its periodic turning points. 
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